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Using turbulent MHD simulations (magnetic Reynolds numbers up to ~ 8000) 
and Hinode observations, we study effects of turbulence on measuring the solar 
magnetic field outside active regions. Firstly, from synthetic Stokes V profiles for 
the Fel lines at 6301 and 6302 A, we show that a peaked probability distribution 



Oh' function (PDF) for observationally-derived field estimates is consistent with a 

O ' monotonic PDF for actual vertical field strengths. Hence, the prevalence of weak 

■ fields is greater than would be naively inferred from observations. Secondly, 
^ ■ we employ the fractal self-similar geometry of the turbulent solar magnetic field 

^ to derive two estimates (numerical and observational) of the true mean vertical 

■ unsigned flux density. We also find observational evidence that the scales of 

(N 

oo 
O 



magnetic structuring in the photosphere extend at least down to an order of 
magnitude smaller than 200 km: the self-similar power-law scaling in the signed 
measure from a Hinode magnetogram ranges (over two decades in length scales 
and including the granulation scale) down to the ~ 200 km resolution limit. From 
the self-similar scaling, we determine a lower bound for the true quiet-Sun mean 
vertical unsigned flux density of ~ 50 G. This is consistent with our numerically- 
based estimates that 80% or more of the vertical unsigned flux should be invisible 
to Stokes— observations at a resolution of 200 km owing to the cancellation of 
signal from opposite magnetic polarities. Our estimates significantly reduce the 
order-of-magnitude discrepancy between Zeeman- and Hanle-based estimates. 

Subject headings: Sun: magnetic fields — turbulence — MHD — techniques: 
polarimetric 



Introduction 



Determining the strength of the magnetization of the "quiet" Sun is tied to the ques- 
tion of how much flux resides at small scales. This is important, for example, in determining 
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the energy budget available for chromospheric heating. Observationally, two methods are 
employed to constrain the solar magnetic field: the Hanle and Zeeman effects. The Hanle 
effect measures (in principle) the mean magnetic field strength, as there are no can- 

cellation effects, but quantitative interpretation requires assumptions about the probability 
distribution function (PDF) of the turbulent magnetic field. The Hanle de-polarization is 
measured in stronger lines formed in the mid- to upper-photosphere an d estimates are made 
o f ~ 130 G, the field residing primarily in the intergranular lanes ( iTrujillo Bueno et al. 



20041 ). The Zeeman effect measures the longitudinal component (via Stokes V) and trans- 
verse components (via Stokes Q and U) of the magnetic field but suffers from cancellation 
effects. Hence, it is "blind" to any "hidden" mixed-polarity flux at scales smaller than the 
resolution limit of an instrument. However, the beneflt of Zeeman measurements is that 
their interpretation requires no assumption about the turbulent PDF and measurements can 
actually be used to determine the PDF on scales larger than the resolution limit. Some 
attempts to incorporate the effects of canc ellation into Zeeman-based St okes inversions uti- 
lize the micro-structured fleld hypothesis (jSanchez Almeida et al.lll996l ). However, typical 
estimates based on the longitudinal Zeeman effect give (l-B^I) ~ 10 G for the mean unsigned 
vertical flux density (see Table 3 in lBello Gonzalez et al.ll2009l for a review of the spread of 
recent Zeeman-based estimates). These values are signiflcantly smaller than the Hanle-based 
estimates. The discrepancy between the results from Hanle (~ 100 G) and Zeeman (~ 10 G) 
measurements is not unexpected since the Zeeman observations see only the resolved flux 
while the Hanle interpretation depends on assumptions about an unknown PDF. 

In this work, we ask the question "How can the turbulent fractal geometry of the mag- 
netic fleld in the solar photosphere be accounted for in interpretations of Zeeman-based 
observations?" We attempt to derive the spectral/fractal properties fr om high resolution 
(0^^3) Zeemari observations with the Hinode spectro-polarimeter (SP) (IKosugi et al.l 120071 : 
Tsuneta et al.l l2008l : iLites et al.l l200ll ) and high resolution (down to 4km;^ 0".006), tur- 
bulent dyna mo sirnulations (up to magneti c Reynolds numbers, Rcm ~ 8000) with the 
MURaM code flVoglerll2003l : IVogler et al.ll2005l ). First, we clarify the consistency of PDFs de- 
rived from observations and simulations, respectively. Zeeman observations typically show 
PDFs of quantities der i ved from Stokes V which can be described as a peaked function 
(IKhomenko et al.l l2005l : iLites et al.l l2008bl ). We demonstrate that the peaked PDFs from 
Stokes V measurements and the monotonic PDFs typically reported from numerical sim- 
ulations are, in fact, compatible. Failure to take into account this observational bias can 
lead to a gross underestimation of the prevalence of weak flelds and, consequently, to incor- 
rect estimates of the mean magnetic fleld strength and magnetic energy density in the solar 
photosphere. 



Second, we extrapolate the results obtained with observations and simulations (based 
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on their respective fractal properties) to scales below the resolution limits to estimate the 
amount of "hidden" flux. To do this, we exploit the fact that the fractal geometry is in- 
timately connected to power-law scaling relations such as N{1) oc l~^f . Here, N{1) is the 
number of boxes of edge length I covering a fractal set (such as all pixels with apparent 
magnetic flux above some threshold) and -D/ is the fractal dimension. It has long been 
known that the distr i bution of plage rn agnetic field is such a statistically self-similar fractal 
( ISchrijver et al.lll992l : lBalke et al.lll993l ). As the threshold can influence the fractal dimension 
inferred, however, the geometrical structure is more complicated than a simple fractal. In this 
case, the fractal concept must be generalized by adding a measure defined by the absolute 
value of the net magnetic flux through each box of edge length /. This measure also displays 
self- similarity (a power-law scaling) for the solar quiet Sun netw ork mag netic flux (down to 
0".5 resolution in Lawrence et al. 1993 and Cadavid et al. 1994: see also Krivova fc Solanki 



20041 ) and also for numerical simulations of magnetoconvection flBrandenburg et al.lll992l ). 
Hence, the geometry of the magnetic field is said to be multifractal. 

Cancellation effects do not play an important role in observations of the multifractal 
unipolar solar regions discussed above. In contrast, if the magnetic field of the quiet Sun 
internetwork is multifra ctal (as we will show), cance llation can play a significant role for 
scales below 1" (see also lSanchez Almeida et al.ll2003l ). For these mixed-polarity fields it is 
necessary to further generalize the fractal con cept to a signed measure. Here, the power-law 
scaling exponent is the cancellation exponent (jOtt et al.lll992l ). If this self-similar power-law 
scaling extends below the observational resolution, small scale cancellation will occur and 
correct values of the rnean fie ld strength and energy cannot be established. As was pointed 
out by [Lawrence et al.l (Il993l ). such a signed measure could also be employed to extrapolate 
moments of the magnetic field below resolvable limits. Until now, no attempt has been made 
to employ self-similar scaling to estimate the total cancellation, and, hence, the true mean 
quiet-Sun magnetic field strength. 



2. Data and Methods 



We use MURaM simulations (IVoglerl l2003l : IVogler et al.l |2004J . l2005l ) for a rectangular 
domain of horizontal extent 4.86 x 4.86 Mm^ and a depth of 1.4 Mm. Runs of small-scale, 
local dynamo action with increasing resolution and Reynolds numbers h ave been carried out 
(see T able [T])- The turbulence is sufficient for small-scale dynamo action (IVogler fc Schiissler 
20071 ). All magnetic field in the simulations results from dynamo amplification of a small 
seed field as there is no net flux through the box, no flux advected into the box, and no 
initial large-scale field. A well-known rule-of-thumb for forced turbulence simulations is that 
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one decade of scales smaller than the forcing is required before the inertial range begins and 
another decade is required (at the opposite end of the spectrum) for the dissipative scales. 
For photospheric magneto-convection simulations, where the forcing granulation scale is 1 
Mm, this means that a grid resolution of 10 krn or srn aller is required before a simulation 
becomes turbulent. (See ISanchez Almeida et al.l (120031 ) for a study with a 10km resolution 
using a Boussinesq simulation and multi-component Milne-Eddington line synthesis.) As we 
compute simulations down to a grid resolution of 4 km, we are now able to measure the 
effects of turbulence on observational quantities. 

To this end, we require synthetic profiles (for the Fel lines at 6301.5 A and 6302.5 A) as 
calculated from a snapshot of a run with non-grey radiative transfer. Owing to the computa- 
tional expense of such a run, it was started from a snapshot from the statistically stationary 
state of Run C (see Tabled]) and then run for approximately one convective turnover time, 10 
minutes. Stokes V, Q, and U profiles were then computed in one dimension (ID) assuming 
l ocal-thermodynamic-equilibrium (LTE) and using the STOPRO code in the SPINOR package 
( 1Solankilll987l : iFrutiger et al.ll2000l ). We will concentrate in this work, however, on Stokes V 
observations (and synthetic observations) at disk center. This allows us to avoid line-of-sight 
effects, to identify the longitudinal component of the magnetic field as its vertical component, 
and to avoid the difficult estimation of the cancellation properties of Stokes Q and U. 



Observational data is obtain ed from the spectro- polarimeter (SP, iLites et al.l 120011) of 



the S olar Optical Telescope (SOT. lTsuneta et al.ll2008l ) on the Hinode satellite flKosugi et al. 



20071 ). One set, a spatial map, consists of 2048 scans taken on March 10, 2007 (11:37:36 
- 14:36:48 UT), in the "normal mode" (exposure time of 4.8s) with the scanning step of 
0.1476" and pixel size along the slit of 0.1585". The data set covers a quiet Sun region at the 
disk center, over the large field of view of 324" x 164". The second set, a "deep mode" time 
series, consists of a 103 steps at disk center, each with an effective exposure time of 67.2 s after 
application of a temporal running mean, and was taken o n February 27, 2007 (00:20 - 02:20 
UT). Both data sets have previously been described in iLites et al.l (l2008bl ). Corrections 
for various instrumental effects are made using the SolarSoft0 procedure sp_prep which 
calculates the wavelength-integrated Stokes V (V^ot); 



j''V{X)dX\ + \J^:V{X)dX\ 



Vtot = sgn(H)- 



Ao 



/ f 



dX 



where sgnfVf,) is t he sign of the blue peak, Aq is the line center, and A,.,;, = Aq ± 30 pm 
(ILites et al.ll2008bl ). The procedure includes also a Milne-Eddington-based calibration of the 



^http: / /www. lmsal.com/solarsoft/ 
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sum of V^r,t from both lines into a measure of longitudinal "apparent flux density", 5^^^ 
( iLites et al.l l2008a[ ). This calibration is tailored to retrieve the field value from weak and 
noisy internetwork signals. It assumes that magnetic structures are spatially resolved (fill 
the resolution element) and does not take into account the magnetic field variations over the 
height range where lines are formed. Consequences of the latter assumption are studied in 
the next section, when -Bfpp, obtained from Stokes V profiles synthesized from simulations, 
is compared with the vertical component of the actual magnetic field. 



3. Results 



3.1. PDFs 



A marked difference exists between the PDFs inferred from Zeeman polarimetry (e.g. 



Khomenko et al 



2005 



Lites et al. 



Socas-Navarro &: Sanchez Almeida 



2003 



2008bl ) and the PDFs frorn nurn erical computations (e.g. 



Vogler fc Schiisslerl 120071 ). For example, in Fig. [T] 



we present the PDFs from the Hinode "normal mode" map magnetogram (apparent vertical 
magnetic flux density, dashed line) and of MURaM simulation Run C-NG (average vertical 
magnetic field, solid line). The polarimetric observation peaks at -Bapp ~ 3G while the 
simulation possesses a monotonic distribution without a distinct maximum: there exists a 
greater amount of weak vertical field than indicated by the observations. The observation 
also shows greater intermittency (the distribution has an enhanced strong signal (field) tail 
when compared to a Gaussian) than the simulation. This can possibly be attributed to 
the much lower Reynolds number of the numerical simulation compared to the Sun as well 
as to the smaller simulation box: contributions from dynamo action in the deeper layers 
and supergranular network flux concentrations are absent. To address the qualitative dif- 
ference (peaked versus monotonic) between simulations and observations, we ask what PDF 
of Zeeman-based observational signatures would result if PDF (5^) monotonically decreases 
with increasing vertical field strength (instead of possessing a peak). Our approach will be 
to assume the distribution of field strengths from turbulent MURaM simulations and examine 
the consequences of such a distribution on Stokes V observations. 

Though noise, resolution, and other instrumental factors are important in any real obser- 
vation, we first address the question assuming a "perfect" instrument. Using the synthetic 
profiles from Run C-NG, we calculate V^q^ with Eq. ([1]) and determine -Bapp following 
Lites et all J2008bh . In Fig. El the derived bI^^^ signal versus -Save? the vertical magnetic 
field strength averaged over the height range corresponding to logr G [— 3.5,.l], is shown. 
This quantity was selected for its linear Pearson correlation with -B^pp of r = 0.92 and its co- 
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efficient of linearity, -Bapp ~ l.O-Bave? wliicli is consistent witli tlie calibration of Lites et al. 



( l2008bl ). This height range also encompasses most of the formation height of the Fel lines 
at 6301 and 6302 A. Though -Bapp and Save are well correlated, there is a large scatter. 
We note, also, that changing the range to logr G [— 2, .1] does not significantly affect the 
correlation, r. This indicat es that most of the Stokes V signal is generated in deeper layers 



( lOrozco Suarez et al.l 120071 ). 

In Fig. [H we present a comparison between the PDFs of -Bapp as derived from the 
synthetic Stokes V profiles (dash-dotted line) and -Bave (solid line). PDF(_Bave) monotoni- 
cally decreases with increasing field strength while PDF(i?app) shows a peak near 1 G and a 
strong decline towards smaller field strengths. The PDFs for maximum Stokes V amplitude 
and total circular polarization are qualitatively similar to that shown for -Bapp. PDFs for 
the vertical magnetic field from different volumes and 2D cross sections from all simulation 
runs, chosen either by height or by optical depth, show similar PDFs to that shown for -Bave- 
That is, the monotonically decreasing distribution is a robust feature of the vertical magnetic 
field when sampled by geometrical height, optical depth, or by averaging over the vertical 
direction. The difference between observations and simulations is caused by the radiative 
transfer that produces circular polarization from longitudinal magnetic field. 

The above result shows that caution is needed when interpreting the distribution of 
Stokes V signal in order to avoid a drastic underestimation of the occurrence of weak field. 
This caution naturally extends to moments of the distribution such as mean vertical fiux 
density or mean vertical magnetic energy density. For example, -Bapp and -Bave are very 

well correlated with a coefficient of linearity of unity, but their averages, ^|-Bapp|^ = 6.9 G 
and (l-Bavel) = 5.5 G, differ significantly. We see that an over estimation of 26% results 
from assuming the vertical magnetic field to have the same distribution as the signal derived 
from Stokes V, even in the a bsence of noise (note that this is close to the 20% loss found in 



Sanchez Almeida et al.ll2003l ). 



To understand in more detail how radiative transfer affects contribute to a peaked PDF, 
we examine a few selected l^— profiles. Pixels with weak -Bave must be generating strong 
-^app signals for PDF(-B|pp) to become peaked. Indeed, Fig. [3] shows this is the case. There 
are many pixels for which |-Bave| < 0.1 G while |-Bapp| > 5 G. On the other hand, we see 
that when |-Bapp| < 0.1 G, |-Bave| is always less than 4 G. In Fig. HJ we examine one case of 
how weak -Bave can be associated with strong -B|p,p. For r G [0.1, 1], the vertical magnetic 
field takes on values of tens of Gauss. In this region, there are also strong gradients (and 
direction reversals) for both the magnetic and velocity fields. In this case, because of the 
magnetic field reversal, -Bave is nearly zero. However, because of the velocity gradient, the 
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contributions to the Stokes— \^ profile from tlie positive and negative magnetic polarities 
are Doppler-shifted with respect to each other. For this reason, the Stokes V signal is not 
cancelled. This is further illustrated in Fig. O where we plot the mean -B^pp for all pixels 
with I -Save I < 0.1 G versus the strength of the vertical velocity fluctuations. There is a clear 
trend of stronger signal with increased Doppler shifts between the different heights in the 
atmosphere. We conclude that Doppler shifts of absorption profiles are responsible for the 
peaked PDF from our noiseless synthetic -Bapp- effect, |V^otl some combination of the 
vertical mean of Bz and the vertical mean of \Bz\ (depending on Vz{z)) and, consequently, 
the PDF of -Bapp does not correspond to FDF{Bz) from any geometrical height, optical 
depth, or volume. Such a failing of -Bapp to accurately represent B^ cannot be captured 

using the Milne-Eddington approximation (used to calibrate -Bapp); which has no gradients 
by definition. 

To examine the effect of noise on the PDF, we consider synthetic Stokes V profiles 
with noise added at a polarization precision of 1.1 x 10"'^ (similar to that of the Hinode 
observations) in determining -Bapp. The PDF of this noisy synthetic observation is shown 
as a dotted line in Fig. [T] and closely resembles the observational PDF for signals weaker 
than a few Gauss. Note that, the noise accentuates the peak in the PDF even further. In 
examining Eq. ([1]) for V^q^ (-Bapp is a nearly-linear function of Vj^ot)' ^^^^ taking 

the absolute value of the blue and red lobes separately the effect of noise becomes the sum 
of two non-negative measurement errors. That is, 

T/measured i/true i i i i /r)\ 

where eb and are the measurement noise in the blue and red lobes (e.g., ef, = JZiti ^i/^ 
where are the random variables associated with the measurement noise in each wavelength 
bin). Assuming these two random variables, eb and have Gaussian distributions, their 
separate PDFs for their absolute values will peak at zero. However, the PDF of the sum of 
their absolute values will peak at a non-zero value due to reduced likelihood that le^l and 
\er\ are small simultaneously: the PDF of the sum of two independent random variables is 
the convolution of their individual PDFs, 

P(6) = r e-(-«)'/2'^^e-«'/2-^rfe (3) 

TTO-feCTr Jo 

for eb,er Gaussian and e = \eb\ + le^]. Because of taking the absolute values, the individual 
PDFs are zero for negative values (this sets the limits of integration for Eq. ([3])). Hence, 
their convolution is zero at zero and peaks instead for some finite positive value. Assuming 



eh and have identical identical standard deviation a = \pl ■ 2.4 G (taken from iLites et al. 



2008b ). it can be shown that the peak in the PDF for -B^pp, Eq. ([3]), is given by the solution 
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to 

B /-^/^ 



a2 



e 



^rfe-e-^'/^ = 0. (4) 



This predicts a peak in the PDF at i?app~ 3 G, close to that the seen in the actual ob- 
servation. Consequently, adding noise leads to a further decrease of the number of pixels 
with very weak field and thus accentuates the maximum of the PDF. This illustrates that a 
monotonic PDF(i?2) is qualitatively consistent with observations. 

In Fig. El we use the 67.2s exposure "deep mode" SP t ime series (dashed line). Note 



that for this exposure time, a = -\/2 ■ 0.6 G (ILites et al.ll2008bl ) and Eq. (jlj) predicts that the 



PDF will peak at -Bapp~ 1 ^ (^^^ actual peak is at ~ 1.2 G). In this case, as in the "normal 

mode" case, the synthetic MURaM -Bapp with equivalent noise level (dotted line) matches the 

location of the peak and the PDF to the left of the peak. In fact, if we generate -Bapp from 
pure white noise for Stokes V (standard deviation of 3 x 10~^), we find its PDF (plus signs) 
predicts well both the location of the peak and the shape of the weak-signal portion of the 
PDFs. This strongly suggests that the observational peak is dominated by noise. 

We also find (see Fig. [6]) that cancellation of opposite polarity fields in a resolution 
element alters the PDF and renders it useless for computations of the mean unsigned flux 
density and other moments. This is evidenced by a comparison of PDFs from the noisy 
synthetic MURaM -Bapp without (dotted line) and wit h spatial degradation by a theoretical 



point spread function (PSF, see iDanilovic et al.ll2008l for details) for Hinode^s optical system 



and rebinned to Hinode pixel size (diamonds). Because of the importance of the cancellation 
on the PDFs, PDFs may not be used to infer the true mean unsigned vertical flux density. 



3.2. Cancellation 

Turbulence gives rise to a statistically self-similar fractal pattern of the magnetic field 
(within the inertial range) - the field retains the same degree o f complexity of distinct 



structures regardless of the sc ale at which it is observed (see, e.g., IConstantin fc Procaccia 



19921 : [Brandenburg et al.lll992l ). In this section, we show that this also applies to solar surface 
magnetic fields and we use this self-similarity to estimate the portion of unsigned vertical flux 
unobservable at a given resolution. To begin with, we examine the cancellation properties 
of the magnetic field itself using a series of high-resolution MURaM dynamo simulations. This 
illustrates how the turbulent nature of the magnetic field limits measurement under the sole 
consideration of spatial resolution and in the absence of other observational constraints. As 
it separates the statistics of the field itself from observational constraints, the study also 
allows us to extrapolate the results to realistic solar magnetic Reynolds numbers. 
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Extending the ideas of sing ularity in probability measures for self-similar fractal fields 
to signed fields, lOtt et al.l (119921 ) introduced the cancellation exponent for studying the self- 
similar sign oscillations on very small scales in turbulent flows. For our application, their 
partition function, xi^)^ measures the portion of the flux remaining after averaging over 
boxes of edge length /, 



xii) 





/a(0 


Bzda 


Ia 


Bz\da 



(5) 



where {Ai{l)} C ^ is a hierarchy of disjoint subsets of size / covering the entire domain, 
A. In our case, we call the function x the cancellation function since it measures the flux 
cancellation at a given length-scale /. If the magnetic field is self-similar (for scales much 
larger than the dissipation scale), we expect a power-law 



X{1) oc l~ 



(6) 



where k is called the cancellation exponent. It is related to the characteristic fractal dimension 
of the magnetic field structures on all scales, Df, by 



{d-Df)/2 



(7) 



wher e ci = 2 is the embedding Euclidean dimension of the solar surface fISorriso-Valvo et al. 



2002|)@ An impr oved method to deter mine x(0 using a Monte Carlo box counting technique 
was proposed by lCadavid et al.l (119941 ) . Its advantages include better counting statistics when 
/ is a large fraction of the edge length of the domain A, applicability to non-square pixels, 
and less sensitivity to the accidental placement of larger flux patches (e.g., network elements) 
with respect to the partitioning. For our simulation data, the Monte Carlo technique proved 
as accurate as rigid partition boxes but led to a significant reduction of the noise in x(0- 
it averages o yer many partitionin gs and allows a more faithful representation of the field 
distribution (jCadavid et al.lll994l ). We use this technique for the results shown below. 



The height range that corresponds (in a horizontally averaged sense) to logr G [—2, 0.1] 
(as discussed in §3.11 the contribution for logr G [—3.5,-2] to the Stokes V signal is in- 
significant) is 2; G [210, 300] km (2 = corresponds to the continuum optical depth r = 1 
at 500 nm). For this height range we compute the averaged cancellation functions, xQ-)-, fo^' 
MURaM dynamo simulations with magnetic Reynolds numbers ranging from Rcm ~ 2000 to 



^Assuming the field is smooth (correlated) in dimensions and uncorrelated in the other d — Df 
dimensions, the smooth dimensions contribute to the sum of vertical fluxes proportional to their area while 
the integral of an unc orrelated field contributes proportional to the square root of its area (random process). 
Eq. d?]) then follows |Sorriso-Valvo et al-lliooi ). 
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Rcm ~ 8000. By definition, we have x(0 = 1 at the resolution of the simulation since there 
are no smaller scales for the computation. Furthermore, we expect dissipation to strongly 
affect x(0 the smallest decade of scales (analytically, its slope must go to zero). Also, as 
our dynamo simulations have zero signed total flux, x(4.86 Mm) = and we would expect 
scales down to approximately 490 km to be affected by this constraint. Only for smaller 
scales should we be able to observe a turbulent scaling. However, very little room is left 
between these two constraints so that no clear power-law scaling is observed for any of the 
simulations (see Fig. [7] for one example). 

Since the dissipation scale of magnetic energy, decreases with increasing Rcm, for 
fixed /, x(0 decreases with increasing magnetic Reynolds number (fluctuations at smaller 
scales increase the total cancellation). This is emphasized in Fig. [HI where we plot the value 
of the cancellation function for / = 200 km (corresponding roughly to Hinode SP's angular 
resolution of 0".3) versus Rcm- We can fit a power law and extrapolate to the results we 
would expect from a M URaM simulation at solar Rcm (which must be estimated). From 



Kovitya &: Cram! (119831 ). we estimate the magnetic diffusivity for logr = 0, 77 ~ lO^cm^s ^. 



The driving of the small-sc ale dynamo i s mainly subsurface where r] is roughly 100 times 



smaller {rj ~ lO^cm^s"^, cf. ISpruitlll974l ). For an upper limit of x(200km), we employ the 
more conservative estimate: 77 ~ lO^cm^s"^. Taking the forcing scale to be the granulation 
scale, £ ~ 1 Mm, and using Vrms ~ 3 km s~^ from the simulation, we find 

Rej,, = ^ 3 . 105. (8) 

For this magnetic Reynolds number, our extrapolation yields x(200km) ~ 0.2. This indi- 
cates that with a perfect observation at this spatial resolution and assuming that the MURaM 
simulation faithfully reproduces the solar conditions, we should multiply an observation by 
a factor of 5 to obtain the true mean vertical unsigned flux density of the quiet-Sun inter- 
network. It is also suggested by Fig. [8] that x(200km) decreases with decreasing magnetic 
Prandtl number, Pm = ^/f] where v is the kinematic viscosity and 77 the magnetic diffusivity. 
As the magnetic Prandtl number of the Sun is much less than that of the simulations, we 
expect that x(200km) < 0.2. 

The cancellation functions for and for i?app inferred from Run C-NG are shown 
in Fig. [71 We see that the two functions are essentially equivalent. This demonstrates an 
excellent correspondence between the cancellation of the field itself and the signal derivable 
from observations (excluding instrumental effects). Therefore, we may take the cancellation 
of -Bapp as a proxy for the cancellation of Bz- This we now do. 

We present the normalized cancellation function, x(0/x(lMm), for the Hinode SP 
observation in Fig. [9] Without knowing the value of the true unsigned vertical flux, the 
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denominator in Eq. ([5]), the value of the cancellation function can only be normalized to 
some arbitrary scale. We find a self-similar power-law over two decades in length scales, 
demonstrating the multifractal geometry of the turbulent quiet-Sun magnetic field. This is 
somewhat surprising as the dominant granulation pattern at scales near 1 Mm might have 
been expected to affect the cancellation scaling. The cancellation exponent of the scaling 
is ft = 0.26 ± 0.01. This exponent predicts a 20% increase in the observed mean unsigned 
vertical flux density with a doubling of resolution in agreem ent with the differ ence in flux 
densities found between ground and space-based telescopes ( iLites et al.ll2008bl ). Note also 
that the power-law behavior holds down to the two-pixel scale. Th is is a clear indication of 
cancellation extending to smaller scales than resolved by Hinode f lCarbone fc Brunol 119971 : 
Sorriso-Valvo et al.ll2004f). Compare this, for e xample, to the simulation case in Fig. [7] (also 



see Fig. 3(b) of ISanchez Almeida et al.l l2003l ) where dissipation is seen to affect a strong 
turnover in %(/) for the smallest decade of scales. As the observation is not so affected, we 
may safely conclude that the smallest scale of magnetic structuring is at least one decade 
smaller than the Hinode SP resolution limit. The scales of magnetic structuring in the 
photosphere must therefore extend to at least an order of magnitude smaller than 200 km. 

From Eq. ([7]), we see that our result corresponds to Df = 1.48 ± 0.02 for the fractal 
dimension of the quiet Sun internetwork magnetogram. W ithin uncertaintie s, this is the 
same dimension as for solar plage regions, Df = 1.54 ± 0.05 (IBalke et al.lll993l ). This might 
indicate that some similar mechanisms are at play in solar plage and quiet S un internetwork 
For t he cancellation ex ponent of network rn agnetic fields, values of k ~ 0.4 flLawrence et al. 
19931 ) and of k ~ 0.12 (jCadavid et al.lll994l ) have been reported, but without an estimate of 
the uncertainties. 

Recent work has highlighted the s ensitivity of fractal d imension (perimeter-area) esti- 
mators to pixelization and resolution (jCriscuoli et al.l 120071 ). By using a signed measure, 
however, we avoid difficulties inherent to fractal dimension estimations using bi-level images 
in general and the perimeter-area method, specifically. Nonetheless, we have tested the sen- 
sitivity of the cancellation exponent to reducing our resolution by theoretical point spread 
functions for apertures 1/2 and 1/4 that of the Hinode SOT (50 cm). We find the slope of 
x{l) to be robust in these cases for lengths exceeding 30 pixels. There is no change in the 
power law for almost one decade of length scales (3-20Mm). We therefore conclude that our 
estimation, k = 0. 26 ± 0.01 is robust and insensitive to pixelization and resolution effects. 
As pointed out by iLawrence et al.l (119961 ) . however, because of what they call "resolution- 
limited asymptotics" , different definitions of fractal dimension can give different values at 
finite resolution. For this reason, our value Df = 1.48±0.02 might differ from a well-resolved 
perimeter- area estimate. 
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Using the self-similar power law derived from Fig. [9l we may estimate the true mean 
unsigned vertical component of the magnetic field (hereafter, "mean unsigned vertical flux 
density") in the quiet-Sun photosphere, (l-B^I). Below the magnetic dissipation scale, 
there is no cancellation: xi^ri) = 1- This, together with the self- similarity relation, Eq. 
gives 

m) = m\ = {m)r(^^J , (9) 

where / is any scale in the inertial range, {\Bz\)i is the mean absolute value of the vertical 
component of the field measured at that resolution (/), 



I A da 



x{i) ■ m) 



and {\Bz\) is given by 



A 



da 



J^da 
11. 7G. As / 



(10) 



'ID 



0.11 Mm (approximate Hinode SP 
L 



Lites et all I^OQM) report (|fi^|)o.ii 

pixel size) is below the SOT resolution limit, however, we rebin -Bapp ^'^ ^ ~ 0.22 Mm pixels 
to find (|i?^|)Q22 ~ 10.7 G as the starting point of our estimate 



(|fi,|)^10.7G 



0.22 Mm 



0.26 



(12) 



Estimating the magnetic dissipation scale is not straight-forward. As we have shown 
that observationally it is unresolved, we are l eft to rely o n a phenomenological estimate. 
Kolmogorov phenomenology predicts (see, e.g., Frisch 1995 ) Ir, ~ CRcj^^^ where £ is a large 
characteristic scale, such as the granulation scale. Using Rcm ~ 3 ■ 10^, derived previously, 
we estimate ~ 80 m. For the dissipative range, power-law scaling for %(/) will not apply 
and the slope of the cancellation function will approach zero. To provide a lower bound 
to the solar mean unsigned vertical field, we should then be conservative by ignoring any 
cancellation in the first decade of scales. Hence, we use = 800 m in Eq. f|T2l) to estimate 
the true mean unsigned vertical flux density to be ^ 46 G. This means that, at a 

resolution of 200 km, at most one quarter of the unsigned vertical flux is observable. 



4. Discussion 



Our estimates suggest that three-quarters or more of the vertical unsigned magnetic 
flux is cancelled at the resolution of Hinode. Hanle-based estimates suggest {\B\) ~ 130 G 
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dTmiiUo Bueno et al.l[20oilN while Zeeman-based estimates suggest (|-Bz|) ~ lOG (see Table 



3 in lBello Gonzalez et al.ll2009l ). Note that even with estimation of the cancellation, there re- 
mains almost a factor of 3 difference between reported Hanle estimates and the Zeeman-based 
estimates we present. However, we have considered only one component of a vector quantity 
while the Hanle-based estima tes are sensitive to the magnitude of that vector. Recent ob - 
servations (ILites et al.ll2008bl ) and simulations (ISteiner et al.ll2008l : ISchiissler &: Voglerll2008l ) 
suggest that horizontal fields are on average a factor of 5 stronger than vertical fields. There- 
fore, our estimate of (l-Bzl) ^ 46 G coupled with an even stronger mean horizontal field is 
consistent with the Hanle-based estimate. Another observational discr epancy lies in deter- 
mining the mean "location" of the fields. iTrujillo Bueno et al.l (120041 ) interpret scattering 
polarization from molecular C2 to indicate the mean field strength is weak (~ 10 G) over the 
bright granules, so that the turbulent field inferred from the Hanle measurements should be 
concentrated in the intergranular lanes. This should be compare d to the location o f strong 
horizontal fields not in the lanes but near the edges of granules (ILites et al.ll2008bl ). In re- 
solving the details over the location and strength of the mean components of the magnetic 
field, future work should also address the cancellation statistics of the horizontal field (and 
the linearly-polarized Stokes signals Q and U) as well as the effect on Stokes V presented 
here. 



5. Conclusion 



On the basis of surface dynamo simulations, we have demonstrated that the PDF gen- 
erated from the Stokes V spectra are not necessarily equivalent in form to that of the PDF 
of the vertical component of the underlying magnetic field. The PDF for Stokes V shows 
a reduction of likelihood for weak vertical magnetic field compared to the PDF of the field 
itself. This effect is not due to a reduction in horizontal resolution, but is caused by a com- 
bination of vertical radiative transfer through a turbulent fluid (via the Doppler effect) and 
noise. That is, any systematic sampling (by geometrical height, optical depth, or volume) of 
Bz from the simulation yields a monotonic PDF, but due to Doppler shifts between different 
atmospheric heights, the Stokes V signal is not such a systematic sampling. Consequently, 
the PDF of Stokes V field estimates do not accurately represent the PDF of the actual verti- 
cal magnetic field even in the absence of noise. Additionally, for two different levels of noise 
( "normal mode" and "deep mode" ) we have demonstrated that the peak in the observational 



■^Rather than as suming that a turbul e nt ma gnetic field possesses a delta-function PDF which leads to the 
^ 60 G estimate in lTruiillo Bueno et all (|2004l ). we take here the ~ 130 G estimate fi-om their assumption of 
an exponential PDF. 
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PDF is dominated by the influence of noise. Because of these two effects, a monotonic PDF 
for the field can result in a peaked PDF in observations and the assumption that PDF{Bz) 
can be uniquely derived from Stokes V observations becomes dubious. 

From the cancellation function for a Hinode observation of the apparent longitudinal 
flux density, we have demonstrated that the multi-fractal self-similar pattern of the quiet- 
Sun photospheric magnetic field covers two decades of length scales down to the resolution 
limit, 200 km. This constitutes observational evidence that the the smallest scale of magnetic 
structuring in the photosphere is at least an order of magnitude smaller than 200 km. The 
power law also allows us to constrain the quiet-Sun true mean unsigned vertical flux density. 
We estimate the lower bound to be ~ 46 G. Estimates based solely on our numerical simula- 
tions suggest that the vertical unsigned flux at Hinode's resolution should be multiplied by 5 
to obtain the true vertical unsigned flux (i.e., ~ 50 G). These two results are consistent and 
suggest that the order of magnitude disparity between Hanle and Zeeman-based estimates 
may be fully resolved by a proper consideration of the cancellation properties of the full 
vector field. 
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Note added in proof 

We would like to po i nt ou t the correlations between our concl usions and the works 
of ISanchez Almeida et al.l (119961 ) and lSanchez Almeida fc Lited f l2000l ) who postulate struc- 
turing of the magnetic and velocity fields on scales much smaller than 100 km. They find 
that synthetic profiles generated by 3-component Milne-Eddington atmospheres re-produce 
the observed Stokes-l^ asymmetries found in 1" resolution observations. Though they did 
not estimate the undetected pho tospheric magnetic flux, th e resul ts indicated a significant 
fraction remaining undetected. iDommguez Cerdena et al.l (120061 ) assumed the quiet-Sun 
PDF can be approximated by a linear combination of the PDF inferred from Zeeman ob- 
servations and a log-normal distribution (accounting for the observed Hanle depolarization). 
They determined that th e Hanle and Zeeman signals a re consistent with a single PDF with 
{\B\) > 100 G (see also ISanchez Almeida et al.ll2003l ). ISanchez Almeidal (120061 ) assuming 
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that numerical simulations of magnetoconvection with no dynamo action (20 km horizon- 
tal resolution) had achieved the asymptotic rate of magnetic energy dissipation, derived 
an estimate for the unsigned magnetic flux contained in unresolved scales; in our nota- 
tion their finding is yf lOOkm) ~ 0.5 while our extrapolation estimates ~ 0.36. Finally, 



Sanchez Almeidal (120081 ) using observational data from various sources plot {\Bz\)i versus / 
(their Fig. 1). The data are c ompared to a lirie, the slope of which corresponds to k = 1, 
i.e., the result for white noise (IVainshtein et al.lll994l . also set Dj = in Eq. ([7])). There is, 
however, a large scatter about this line suggestive of either a large uncertainty in n or an 
element of randomness in the calibration issues between the various data used. 
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10"^ ^qO ^q2 ^q3 

B (G) 

Fig. 1. — Probability distribution functions (PDFs) for magnetic field strengths and derived 
field proxies: Hinode SP "normal mode" map -Bapp (dashed line), MURelM simulation -Bave 
(see text, solid line), MURaiM synthetic -Bapp (^ derived from Stokes V, dot-dashed), and 
-^app including a noise level of 1.1 x 10~^ (dotted). The PDFs of the synthetic observations 
appear peaked although we have a monotonic distribution of vertical field strengths. 
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Fig. 2.— B^pp derived from MURaM Run C-NG versus -Save? the actual vertical magnetic 
field strength averaged over logr G [— 3.5,.l]. The linear Pearson correlation for the two 
quantities is r = 0.92. Note the large scatter. 
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Fig. 3.— (Left) ^Lpp versus B^ve for B^ve < 0.1 G (Right) B^ve versus E^pp for 5^ 



< 0.1 G (-Bapp computed from noiseless V-profiles). These pfots indicate the bias that 
strong Stokes V signal can be associated with a pixel with weak averaged magnetic field, 
but seldomly vice-versa. 



app 



'app 




i-og ^ X 

Fig. 4. — (Left) Bz (solid line) and Vz (dashed line) versus optical depth, rsoonm, and 
(Right) Stokes V profile for the pixel indicated by a diamond in Fig. [3] (-Bapp = —5.9 G 
and -Bave = —1-6 ■ 10^^ G). At log r = the positive and negative contributions to -Bave have 
nearly cancelled (integrating downward). The Stokes V signal is stronger than would result 
from a uniform 1.6 ■ 10~^ G field but is asymmetric. Strong gradients lead to asymmetric 
profiles but also to I-Bappl ^ I -Bave |- 
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Fig. 5. — Average i?app versus standard deviation of the fluctuations of the vertical velocity 
along the (vertical) line-of-sight, c^^,, for all pixels with |i?ave| < 0.1 G. Pixels are binned by 
(Ty before averaging. With strong velocity differences between different heights in the atmo- 
sphere, the total Stokes V signal increases as the Doppler-shifted absorption from positively 
and negatively oriented fields show less cancellation. 
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Fig. 6. — PDFs for derived field proxies: Hinode SP "deep mode" time series -Bapp (dashed 
line) and MURaM synthetic -Bapp including a noise level of 3 x 10~^ (dotted). The effects 
of cancellation due to finite spatial resolution are seen in the PDF of the synthetic signal 
including noise as well as spatial smearing from a theoretical PSF and rebinning to Hinode 
resolution (diamonds). As this represents a real loss of data, the true mean unsigned vertical 
flux density cannot be calculated from the observational PDF. Also shown is the result for 
employing pure white noise with a standard deviation of 3 x 10~^ for Stokes V (+). 
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Fig. 7. — Cancellation function, xi^)^ versus scale, / for Run C-NG: (solid line) and 
-^app (simulated observation, dashed). The two are essentially equivalent, suggesting that 
the cancellation of -Bapp may be taken as a proxy for the cancellation of Bz- 
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Fig. 8. — Portion of flux remaining at / = 200 km, x(200km), versus magnetic Reynolds 
number, Rcm- Symbols are Run E (plus). Run C (asterisk). Run G (diamond). Run H 
(X), Run C-NG (triangle), and Run G-P (square)-see Tabled! For fixed /, x(/) decreases 
with Rcm and shows an approximate power-law relation with Rcm as indicated by the fitted 
dashed line. Run C-NG and Run G-P are not included in the fit, but the effect of decreased 
magnetic Prandtl number leads to reduced x(200km). Taking this into consideration, along 
with extrapolation to solar values, Rcm ~ 3 ■ 10^, we estimate x(200A;m) < 0.2. 



- 25 - 




0.1 1.0 10.0 

I (Mm) 

Fig. 9. — Normalized cancellation function, Mm), versus scale, /, from Hinode -Bapp 

observation. A self-similar power-law is abundantly clear for 2 decades of length scales down 
to the resolution limit of the observation (the fitted line is k = 0.26 ± 0.01). This indicates 
both the possibility for self-similar extrapolation to smaller scales and that the smallest 
scales of magnetic structuring must be at least an order of magnitude smaller than 200 km. 
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Table 1. Summary of MURaM simulation runs: shown are grid points, horizontal resolution, 
and magnetic Reynolds number, Rcm- All runs exc ept Run C-NG util ize grey radiative 
transfer. In Run C-NG, opacity binning with 4 bins (jVogler et al.ll2004l ) has been used to 
provide non-grey radiative transfer. For all simulations no physical viscosity is imposed. 
R ather, numerical di ssipative effects lead to an effective kinetic Reynolds number. Re 
( Vogler et al.lboosi ). To obtain a lower value of Pm = Rcm/ R^, Run G-P uses the 
magnetic diffusivity used in Run C but at a higher resolution, hence higher Re. 
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